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The effect of a small amount of noise on the standard mapping is considered. Whenever 
the standard mapping possesses accelerator modes (where the action increases approxi- 
mately linearly with time) , the diffusion coefficient contains a term proportional to the 
reciprocal of the variance of the noise term. At large values of the stochasticity param- 
eter, the accelerator modes exhibit a universal behavior. As a result the dependence 
of the diffusion coefficient on the stochasticity parameter also shows some universal 
behavior. 
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1. 



Introduction 



The standard mapping, 



Vt - Vt-i = -esmxt-i 



Xt - Xt-l = Vt, 



(1) 



is an important model with which to study the phenomenon of stochasticity [1] . Recently, 
Rechester and White [2] calculated the diffusion coefficient for this mapping with an added 
noise term Sxt in the equation for xt — Xt-i- The random variable Sxt is sampled out of a 
normal distribution f{dxt;a) of variance a, 



Noise has been introduced into the standard mapping by other authors [3, 4]. However, 
they have primarily considered e < 1 where the standard mapping (without noise) has large 
regular regions. The importance of the noise is to allow a trajectory to wander from the 
regular to the stochastic regions leading to a large change in the quality of the motion even 
when the noise is small. 

We shall be concerned here with the effect of a small amount of noise on (1) in the 
stochastic regime e > 1. It is not so clear in this case how important noise will be. The 
primary effect of the noise is to make the motion ergodic. "Ergodic" here means that for 
a single realization of the noise, time and phase-space averages will be equivalent. (When 
taking averages, phase space is defined by taking x and v modulo 27r.) In the stochastic 
regime, there are often small regular regions present with a = 0. Since their total area is 
small, one might think that they have little effect when noise is introduced. However, if some 
of the regular regions are "accelerator modes," which are stable regions in which the particles 
are continually accelerated [1], the diffusion coefficient may exhibit a 1/cr dependence for 
small a. We may easily see the origin of this 1/a dependence. Ergodicity ensures that a 
fixed fraction of time is spent in each accelerator mode. As a is decreased, the mode is 
visited less often while the duration of each visit is proportionately increased. (Only the 
noise can cause a particle to enter or leave an accelerator mode.) This leads to a diffusion 
consisting of a few large steps. Since the step size Av and the duration of the step are both 
proportional to I/ct, the contribution to the diffusion coefficient also scales as l/cr. In this 
paper, we shall confirm the 1/cr dependence for two different types of noise. However, we 
expect similar behavior for any noise model. 

The remainder of this paper is organized as follows: In section 2, the results of rcf. 2 arc 
interpreted in terms of the correlation function. This allows us to extend the result of ref. 2 
to include other noise models and to determine the behavior of the diffusion coefficient when 
(7 is small. In particular, we predict the 1/cr dependence when there are accelerator modes 
present. Section 3 looks at the behavior of the diffusion coefficient numerically. The 1/cr 
dependence is confirmed and the dependence on e is found. In section 4, we examine the e 
dependence in more detail. We find that the accelerator modes, which give the 1/cr behavior, 
have a universal structure when e is large. This means that the diffusion coefficient exhibits 
a dependence on e which has some universal characteristics. The results are summarized in 
section 5. 



/(a;;a) = (27r<7)-V2exp(-i 
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2. The correlation function 

We shall examine the standard mapping (1) with two different types of noise. With the 
first type, the standard mapping becomes 

vt - vt-i = -tsmxt-i + 5vt-i, Xt-Xt-i=Vt + Sxt, (2) 

where Svt and 5xt are random variables sampled from distributions f{Svt; p) and f{5xt;cr) 
respectively. This is a simple generalization of the system treated in ref 2 so that the noise 
causes diffusion in both the v and x directions. 

The effect of the second type of noise is described by 

vt- vt-i = -esinxt-i, zt-zt-i=Szt, Xt - Xt-i = Vt - Zt + 6xt. (3) 

Here {Sxt, 5zt) is a random point chosen with a distribution g{dxt, Szt; a) where 

g{x, z; a) = exp(— cr) S{x, z) + [1 — cxp(— cr)]/ (27r)^ u{x, z), 

(5(.T, z) is the Dirac delta function, and u{x, z) is a uniform function equal to 1 if \x\ < tt 
and |z| < TT and otherwise. This noise term models the effect of large-angle scatterings. 
The effect of z in (3) is to provide an origin shift to v in the equation for xt — xt-i- The 
solution to (3) behaves as that of the noiseless standard mapping (1) as long as (Sxi.Szi) 
is (0, 0). After an average of 1/a (for small a) iterations, a large-angle collision takes place 
which completely randomizes the particle's position in x and z. 

Equations (2) and (3) represent examples of two distinct types of noise. In (2) the noise 
is diffusive; this is the limit in which the particle suffers frequent but small uncorrelated 
kicks. On the other hand, in the large-angle scattering model (3), the particle is rarely 
kicked by the noise but the kicks are large. Chirikov [4] has considered (2) when cr = 0, 
although he introduced the noise as a shift of the origin of v rather than directly into the 
V equation. [Thus Chirikov's model is given by (3) with Sxt = and dzt sampled from 
fiSzup).] 

In the absence of noise, there exist small regular regions from which an orbit is excluded 
(assuming it started outside such a region). Noise destroys these regular regions allowing 
a particle to wander anywhere in phase space. In addition, the types of noise we consider 
in (2) and (3) maintain the area-preserving nature of the standard mapping because at any 
given time step they merely translate the phase space by some fixed amount. Prom this it 
follows that the motion is ergodic, i.e., that time averages can be replaced by phase-space 
averages. 

The velocity-space diffusion coefficient is defined by 

t— >oo zt 

where the angle brackets denote an average over some appropriately chosen ensemble. An 
equivalent definition is [7] 

oo 

D = i^Co + Y,Cr, (5) 

where 

at = Vt+1 — Vt is the acceleration, and the average now includes an average over t. Because 
the motion is ergodic, we can replace the time average by a phase-space average coupled 
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with an average over all realizations of the noise terms. (In calciilating phase-spac;c; averages 
we use the periodicity of the mappings in the x and v directions so that the averaging need 
only be done over a 27r x 27r square. However, when defining D and a^, the periodicity in v 
is not used.) The result for will be independent of the ensemble chosen so the ensemble 
average can be ignored. Thus we have 

Cr = ( / TT^ / TT^arOo ) , (6) 



/o 27r Jo 2-K 

where here the angle brackets mean an average over the distributions of all the noise terms 
appearing in aT-aQ. For instance, for the noise term in (2) we have 



/oo /"CO 
-oo J —oo 

/oo pOO 
f{Sxi;a)d5xi / f{dvo;p)d6vo 
-oo J —oo 



This operation is an identity in the limit a and p — » 0. 
The first few Cr are then found for (2) to be 

Co = \<? + P, Ci = 0, 

Ci = -ieV2(e)exp(-CT- \p), 

Cs = -\(^Jl{€) cxp(-CT - p) + ieV^e) exp(-3CT - p). 

Assuming that the sum in (5) can be tnmcated at t = 3, we have 

D w ip+ ie2[i - J2(e)exp(-(T - \p) - (e) exp(-cr - p) + j|(e) exp(-3cr - p)]. 

With p ^ this agrees with the result obtained by Rechester and White [2]. Similarly 
when cr 0, D — ^p is the diffusion coefficient given by Cohen and Rowlands [5] for 
Chirikov's noise model [4]. (We must subtract ^p to account for the different way in which 
noise is introduced by Chirikov.) We note that the expression in Ref. 5 also includes a 
term ^e^ JKe) exp(— p) which is part of C4. We might expect the truncation of (5) to be 
accurate when e greatly exceeds the stochasticity threshold, i.e., e ^ 1. This question will 
be examined in more detail below. This approach shows that the oscillations in D seen by 
Chirikov [1] are due to short-term correlations in the standard mapping and, contrary to 
his assertion, are not directly caused by the presence of accelerator modes. 
The same calculation may be made for (3). Here we obtain 

Co = , C\ = 0, 

C2 = -ieV2(e)exp(-2a), 

C3 = H^^Jli.^) + ieVf(e)]exp(-3<7), 

Cr = exp(— crr)CT-(<T 0). 

This result for Ct comes about because the probability that at least one large-angle scat- 
tering takes place between oq and is 1 — exp(— ctt). If no scattering takes place then the 
mapping is the same as the noiseless one. 

A similar method for deriving the diffusion coefficient is given by Cary et al. [8] who 
define a more general correlation function. While this method is more complicated than that 
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described above, it docs allow the calculation of other statistical properties of the mapping. 
In order to illustrate this method of computing the diffusion coefficient for another mapping, 
we consider the mapping obtained for the motion of an ion in a lower hybrid wave [7]. The 
mapping is 

xt - xt-i = 2TrS - 2'!rAcosyt-i, yt - Vt-i = 2tt5 + 2TTAcosxt 

and we are interested in diffusion in the v direction where v = ^{y — x). We take A to be 
much larger than the stochasticity threshold A j. Then the motion is approximately 
ergodic (even though there is no noise in this model). Furthermore, we expect (subject to 
the restrictions to be explored later in this section) that only a few terms in (5) contribute 
to D. So we have 

D « tt'^A^I^ + Joi2nA) cos(27r(5) - J^{2ttA) sm^{2TTS)], 

where the first term in the brackets is the contribution from Co and the other terms come 
from Ci. Antonsen and Ott [9] have also derived this result using the method of paths in 
Fourier space [6]. 

The question of the accuracy of discarding the terms for r ^ 3 in (5) may most easily 
be addressed with the noise model employed in (3) because we need only determine the 
behavior of Ci-(a — > 0) for r > 3. When a = and e > 1, phase space may be divided 
into two regions: a large connected stochastic region and those parts of phase space within 
islands. A particle starting in either region stays forever in that region. It is useful to write 
Cr as the sum of C^* and C** which are the contributions to the integrals in (6) due to the 
stochastic and island regions respectively. 

We shall assume that C^* decays exponentially with increasing r. So far as we know, this 
has not been proved for the standard mapping. The numerical evidence is that C** decays 
quite rapidly for small t and large e. The decay for larger r is difficult to measure accurately 
because the error in the measurements of C^* may exceed the value of Cf. (Grebogi and 
Kaufman [10] have numerical evidence confirming the exponential dependence for quite large 
T. On the other hand, the work of Channon and Lebowitz [11] on the Henon map [12] shows 
an algebraic decay of correlations for a stochastic orbit. But their results do not preclude 
an exponential decay for longer times than they were able to measure.) 

The contributions due to the island region may be evaluated quite accurately because, 
within a given island, the time-averaged acceleration is a constant {i is a subscript labelling 
the various islands). The frequency of oscillation around an island is typically of order unity. 
Therefore C** consists of a mean part (independent of r) equal to Q = afAi/Ao plus a 
part which oscillates with a frequency of about unity. Ai is the area of the ith island and 
^0 = (27r)^ is the total area of phase space. 

We are now in a position to assess the contributions to D of C** and C*" when a is 
finite. The sum over C** probably converges rapidly so that a truncation at some fairly 
low T is quite accurate. Since, in that case, only the terms for small r contribute to D, a 
small amount of noise has little effect on this contribution. The oscillatory part of Of may 
similarly be neglected when evaluating D since its sum when weighted by exp(— ar) is on 
the order of a. The mean part of C^'^ , on the other hand, increases D by 
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= 2 + E ^ Q ( 2 + E exp(-aT) 




= iQcoth(ia) w g/cr 




(7) 



This will be nonzero if at least one of the islands is an accelerator mode, i.e., ^ for some 
i. In these modes, a particle, instead of returning to the original island after N iterations, 
goes to the image of that island displaced upwards or downwards in v by some multiple of 
2n. Such a mode is called an A^th-order accelerator mode. (In fact, due to the symmetries 
of the standard mapping, they come in pairs with = ^2i-i and a2i = —a2i-i-) So, if 
accelerator modes exist, a can be chosen so that Dig and hence D are arbitrarily large. 

Accelerator modes are best found by looking for stable accelerating fixed points. Around 
each such fixed point there will be an accelerator mode. Several first- and second-order 
accelerating fixed points for (1) are cataloged in table I. Chirikov [1] gives 



with n being an integer as the condition for the stability of first-order fixed points. The 
magnitude of the acceleration of the accelerator mode associated with such fixed points is 
27rn. When e = ei, the first order fixed points become hyperbolic with reflection and a 
pair of stable second-order fixed points (marked by asterisks in the table) are born. The 
accelerating regions around these second-order fixed points are best thought of as being a 
continuation of the first-order accelerator mode. 

The physical explanation of the 1/cr divergence was given in section 1. Here we will cast 

that explanation into more quantitative terms. Consider a particle that has just been placed 
in an accelerator mode by a collision. When a is small, the probability that it survives in 
that mode for longer than a time t is P{t) = exp(— at). (This is just the probability that 
there is no collision during the time t.) The probability that it leaves between times t and 
t + dt is p(t)dt where p{t) = —dP{t)/dt = C7exp(— at). From (4), the contribution to the 
diffusion coefficient from the accelerator modes is 



The factor Ai/Ao is the fraction of time a particle spends in the ith accelerator mode. 
Substituting for p{t), we have Dig = Q/a which agrees with (7). 

The same considerations apply to the standard mapping with noise given by (2). In 
this case, the form of p{t) is not known; it will in fact depend on the size of the accelerator 
mode. However, we do expect the duration of a particle's stay in an accelerator mode to 
be approximately A^/cr (for p = 0) where A is the scale length of the island. This is to be 
compared with an average duration of 1/a for the large- angle scattering case (3). Since A 
is usually quite small, the coefficient of the 1/a term in D for (2) should be smaller than 
that for (3). These considerations will be refined in section 4, when we will be able to make 
more accurate scaling arguments. 

In the next section we numerically confirm the 1 /a dependence and explore the depen- 
dence of I? on e. 



eo = 2nn < \e\ < [{2nnf + 16]^/^ 



= ei 



(8) 




(9) 
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3. Numerical evaluation of the diffusion coefficient 

In order to measure the diffusion coefRcient numerically wc adopted a method based 
on (5) which is designed to handle systems with long correlations. The trajectories of J 
particles with random initial conditions (with a uniform distribution) are advanced iot = kT 
according to cither (2) or (3). A correlation function is defined for the fcth iterate of the 
map for each trajectory by 

T-r-l 

where a^ = vi.{tj^i) — Vkt is the acceleration due to k iterates of the map. A diffusion 
coefficient based on the jth trajectory is given by 

r=l 

The final value of D is obtained by averaging Dj over j. The standard deviation of Dj 
divided by y/j is used to give a measure of the error in D. 

In this method, correlations up to a time separation of kL are retained. With T = 1 
and i = 0, we recover the "standard" method which is based on (4); in order to obtain 
accurate results in this case J vaxist be large. Here we do not take J to be large; however, 
good statistics are obtained because we take T ~ 1 so that there are many observations 
of each . Normally we take J to be 64 which allows us to make full use of the vectorization 
capabilities of the Cray-1 computer on which the computations of D are performed. 

In fig. 1, we show the a dependence of D for (2) with p = and e = 6.6 and 12.8. D is 
normalized to its quasi-linear value Dqi — ^e^. These values of e were chosen to satisfy (8) 
for n = 1 and 2. We see that D does have a l/cr dependence for small a. The values of a 
at which this dependence becomes evident are about 10"* and 10~^ for e = 6.6 and 12.8. 
The values for D iov a = 10~^ greatly exceed the numerical values given in ref. 2. Since 
orbits of length 50 were used in those computations, the effect of the accelerator modes was 
largely suppressed. 

Also shown in fig. 1 is D for (3) with e = 6.6. As expected, the coefficient of the l/a 
term is nearly 100 times larger than for (2). 

Taking the limit D{<t 0) gives an infinite result from the I/ct term. If we interchange 
the limits so that we take a ^ before t ^ oo in (4), the value of D depends on how the 
initial conditions are chosen since the motion in this case is not ergodic. If an ensemble is 
defined by choosing initial conditions uniformly in phase space, D is infinite because some 
trajectories will be accelerating. This is then consistent with the value of D obtained by 
taking the limit cr ^ after t ^ oo. A more "natural" ensemble is obtained if we restrict 
the initial conditions to the stochastic region of phase space. Figure 1 shows the value of 
the D for e = 6.6 and (7 = with such initial conditions. The error in this measurement of 
D is quite large even though long trajectories were used in the computation. This probably 
arises because the stochastic region includes a "sticky" portion close to the accelerator 
modes. A particle which wanders into this portion of phase space may still spend a long 
time accelerating even though the trajectory is still stochastic. The properties of these sticky 
regions around islands need more thorough study if the diffusion coefficient for a = is to 
be understood. 
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No enhancement of D was detected in the ranges of e where; second-order accelerator 
modes exist. Because these modes are much smaher than the first-order accelerator modes 
the value of a at which they begin to contribute significantly is so small that prohibitively 
long runs would have to be made to detect any effect numerically. 

Next we turn to the behavior of D as a function of e. Here we hold cr fixed and equal 
to 3 X 10~^ and e is varied in and somewhat beyond the ranges given by (8) with n = 1 
and 2. The results are shown in fig. 2. D rises quite rapidly as soon as e exceeds eo for the 
first-order fixed point. At about one quarter and at about one half of the way through the 
interval (eo, ei), D is dramatically reduced. As we shall see this is due to the appearance of 
fourth- and third-order resonanc;es. Nothing much happens to D at ei. Although the central 
fixed point becomes unstable at this value of e, there is still a KAM surface of the original 
topology surrounding both the unstable first-order fixed point and the two new second-order 
stable fixed points. There is little change in the overall size of the island at this transition. 

Perhaps the most noteworthy feature of fig. 2 is that the plots for both n = 1 and 2 are 
so similar. This suggests that there may be a universal structure for the accelerator modes. 
We pursue this subject further in the next section. 

4. Universal behavior of accelerator modes 

Referring to table I, we see that the accelerator modes exist only in a narrow range 
in e. They are likewise present only in a small region of phase space. This effect becomes 
more pronounced as e is increased and allows us to approximate the accelerator modes by 
a Taylor-series expansion of the mapping. We consider a general area-preserving map of 
the {x,y) plane which depends on a parameter k. We shift the origin and k so that the 
accelerator mode first appears at {x,y) = and k = 0. We pick a frame traveling with 
the acceleration of the mode; therefore the constant terms which represent the acceleration 
are subtracted. We shall only directly treat accelerator modes which appear as a result of 
tangent bifurcations. Other higher-order fixed points which come from bifurcations of an 
existing accelerator mode will be treated as part of that accelerator mode. The linear terms 
of the mapping at fc = have the form 

xi = (1 + a')xo + byo, yi = b'xo + [1 - a')2/o, 

with a'^ -|- bb' = 0. Because the mode appears as a tangent bifurcation, the trace of the 
tangent mapping matrix is 2. By transforming (a;, y) with 

x = x' + a'y'/{a'^ + 6^), y = -a'x'/b + by'/{a'^ + 6^), 

the linear mapping becomes 

x[ =x'o + y'o, y[ = y'^. 

A similar transformation is possible if 6 = but b' ^ 0. The only case where the transfor- 
mation is not possible is if a' = 6 = 6' = in which case the linear mapping is an identity. 
We now add the terms in the Taylor expansion which are quadratic in x and y and linear 
in k, 

A = x'a + l/a + ('-''a + d.x'„yo + ey'^ + fk (10a) 
y[=y'o+ dx'S + d'xy^ + e'y';} + fk. (10 b) 
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Wc have taken k ^ O(x^). Unfortunately, this mapping is not in general area-prc;serving 
because the omitted cubic terms in the map also contribute to area preservation. However, 
we can make (10) conserve area by expressing four of the coefficients of the quadratic terms 
in terms of the other two. The coefficients of x'q, c and c', are taken as the independent 
ones. This choice is motivated by noting that the results are then independent of our choice 
of the direction of time. Thus when we invert (10), the coefficients of xf depend only on c 
and c', whereas the other coefficients are linked in a more complicated way. Also, in some 
respects we may order y' as a;'^ [e.g., consider the positions of the fixed points of (10)]. The 
terms involving x'qIJq and ijq are of the same order as the neglected terms. Therefore their 
coefficients should only be chosen to ensure the preservation of area; i.e., 

d' = 2(c' - c), e' = (c' - cf/c', d/d' = e/e' = c/c'. 

Equation (10) may then be written as 

x[ =x'o + y'o + c[c'x' + (c' - c)yf/c'^ + fk, 

y[=y', + [c'x' + {c'-c)y'f/c' + f'k. 

Finally, we perform the transformations 

x' = 2[x" - (1 - c/c')y"]/c', y' = 2y"/c' 

and 

x"=X, y" = Y + 2K{f/f' -c/c^), k = -4K/{c'f) 

to give 

Yi-Yo = 2{X^-K), Xi-Xo = Yr. (11) 

This is a universal mapping approximating the behavior of accelerating modes for large 
stochasticity parameter. All the transformations which give (11) are linear, and, with the 
exception of the last one which is just a shift of the origin, they are all independent of k 
(and K). For the first-order accelerator modes for the standard mapping (1) which appear 
at e = 27m, the transformations reduce to 

x = ±(- + —x\ v = ±(-2Trnt+—Y], e = 2Trn+—K. (12) 
\2 7rn y \ im J irn 

The transformations leading to (11) are not well defined if the various coefficients satisfy 
unusual relationships. The first transformation is not possible if the linear term is an identity. 
If the motion around a stable A'^th-order fixed point of a mapping is similar to rotation by 
an angle 2'Kp/q where p and q are integers, then the linear part of the qNiYi iterate of the 
map is an identity. But this does not correspond to the first appearance of an accelerator 
mode. (It appeared with the A'th-order fixed point or sooner.) The transformations also 
fail if c' or /' is zero. In this case higher order terms have to be kept. This is what happens 
with the standard mapping at its (non-accelerating) fixed point which appears at e = and 
{x,v) = (0,0). The same thing happens if we look at the second-order fixed points which 
appear when a first-order fixed point goes unstable (the fixed points marked with asterisks 
in table I). But here again such second-order fixed points are not the first occurrence of an 
accelerator mode in the neighborhood of parameter and phase space. They are treated by 
the second-order fixed points of (11) that appear for > 1. Another form of degeneracy 
occurs if the exact mapping has, for instance, a real square root. Then quadratic terms in 
the mapping may be finite but still the cubic and quartic terms may not be neglected. An 
example of this is provided by the second iterate of (11). 
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Figurc 3 illustrates this mapping for a particular value of K and also shows that it 
does indeed closely approximate the standard mapping near first- and second-order accel- 
erator modes. Thus accelerator modes may be studied by examining (11) and how the 
transformations affect the diffusion. 

We begin by cataloging some properties of (11). For Q < K < 1, (11) may be trans- 
formed into the Henon quadratic map [12], 

Xx = xq cos a — {uq — x^) sin a, yi = xq sin a -|- [yo — x^) cos a. 

K is related to Henon's parameter K = sin^(a/2). The transformation between the two 
sets of coordinates is given by 

X = cos2(a/2) sin(a/2)«; -|- Y/2 - sm'^{a/2), Y = -2 cos(a/2) sm^{a/2)v, 

and 

w = xcos{a/2) + ?ysin(a/2), v = — a;sin(a/2) + ycos(a/2). 

The transformation depends intimately on a (and hence on K) and becomes singular at 
ii' = and 1. 

For < 0, (11) has no fixed points. It has two first-order fixed points for > 0, 
(X.Y) = (=F\/X, 0), which following ref. 12 we call Ii and /( respectively. I{ is always 
hyperbolic and so is unstable. Ji is elliptic and therefore stable for K < 1. At K = 1, Ii 
turns into a hyperbolic point with reflection and two second-order fixed points are born. 
Thus the values K = and K = 1 correspond to the parameter values eo and ei listed for 
the fixed points without asterisks in table I. 

The second-order fixed points are stable for 1 < < | . This corresponds to the range 
in e for which the fixed points labelled by asterisks in table I are stable. At K ^ j, a. second 
bifurcation takes place giving rise to periodic trajectory of period 4. This trajectory in turn 
becomes unstable at K = 1.2801 when a period-8 cycle is born. The process of a period-2'" 
trajectory becoming unstable and producing a period-2'"+^ trajectory continues. Greene 
et al. [13] show that it accumulates at K = 1.2840 = i(l + 1.2663)^. When K exceeds 
this value, they conjecture that (11) has no stable fixed points. Therefore, for large n the 
first-order accelerator modes exist for 

eo = 27rn < e < £2 « eo 1.2840(ei - eo) w 27rn + 1.2840 x 4/(7rn). (13) 

The map (11) has two main symmetry lines. They arc the X-axis {Y = 0) which 
corresponds to the w axis in the Henon map and the line Y = K — X"^. Reflection in one of 
the symmetry lines corresponds to reversing time. E.g., F = is the invariant line for the 
transformation X = X' — Y', Y = —Y' which turns (11) into its inverse. 

In order to understand how (11) contributes to diffusion we must determine how large a 
region is trapped around the stable fixed points. We define an orbit to be trapped if X and Y 
remain bounded for all time. (Trapped orbits are the ones which contribute strongly to the 
diffusion because in the original mapping they are the ones that are perpetually accelerated.) 
Untrapped orbits escape to infinity with Y +00. It is straightforward to show that for 
K < 2 all particles with F > 10 escape in this way. An approximate numerical test for 
being trapped is to check that y < 10 during a large number of iterations. Figure 4 shows 
the extent of the trapped region along the two symmetry lines and its total area as functions 
of K. Figure 4b is just a rescaling of Henon's fig. 2. The dips in the plots of the diffusion 
coefficient (fig. 2) correspond to the large reductions in the size of the trapped region that 
occur for certain values of K. These are associated with the occurrence of higher-order 
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rcsoiianccs to the basic rotation about the Ii [12]. Particularly large effects are produced by 
the fourth- and third-order resonances. A stable period-four cycle exists for ^ < /sT < 0.3044. 
This is born at Ii and in the process of moving away from this point destroys much of the 
stable region. A stable period-three cycle /a exists ioi ^ < K < — 0.5625. Unlike 
other resonances this is born away from the associated lower-order fixed point /i . moves 
further away from I\ as K increases, while its unstable twin /g moves towards I\. When 
Is goes unstable at JsT = /g is at Ji causing the destruction of all stable regions in the 
neighborhood of Ji. These two resonances cause the large reductions in D seen in fig. 2. 
Thus, the values of e at which the dips occur are given hy K 0.3044 and 0.5625 (with 
n = 1 and 2) in (12). 

In order to make a more quantitative comparison between the behavior of the diffusion 
coefficient in fig. 2 and the behavior of (11), we must ascertain the effect of the transforma- 
tions used to derive (11). We saw in section 2 that the contribution of an accelerator mode 
to the diffusion coefficient when the noise is due to large-angle scattering (3) is given by the 
area of the mode and its acceleration. Now the area of the accelerator mode is 4/(|6| c'^) of 
the area of the trapped region of (11) (fig. 4c). So in the limit of small a, the diffuwion coef- 
ficient for (3) consists of a superposition over accelerator modes of forms of fig. 4c, linearly 
scaled by appropriate amounts in both directions. For a first-order accelerator mode, Ai 
the area of the mode scales as 1/n? ~ 1/e^ (12) while the acceleration varies as Oj ~ n ~ e. 
Thus from (7) we have Dis ~ 1/cr. The relative importance of the first-order modes with 
this noise model is given by Dig/Dqi ~ l/(cre^). 

With the noise employed in (2), the situation is more complicated because the way 
diffusion due to a and p enters (11) depends in a more involved way on the transformations 
used to derive (11). In particular, although fig. 4c has the same quaUtative features as 
fig. 2, the diffusion coefficient is no longer simply a function of the area of the accelerator 
mode. For instance, take the case where /? = (then all the diffusion due to the noise 
is in the x-direction). In (11) the diffusion is then along a line Y = a'{a'^ + h'^)X/b'^. 
So the effect of each accelerator mode as a function of K depends on what this direction 
is. However, in (2) the important accelerator modes are the first-order ones and for these 
modes the transformation to [X^Y) takes the particularly simple form (12). The only way 
the transformation changes as n is changed is by an overall scale factor. The effect of the 
noise will then only depend on the ratio of /? to cr (apart from a constant factor). 

In order to find the effect of a first-order accelerator mode in (2) on the diffusion 
coefficient, we go back to the expression for D^g given in (9). Now p{t) is no longer a 
probability since it can not be normalized. It may be defined as follows. Imagine starting 
infinitely many particles at the boundary of the trapped region. (The number of particles 
has to be infinite because almost all of them leave the trapped region immediately.) Suppose 
the number which are left after a time t is P{t). The number which leave between t and 
t + dt is th.eiip{t)dt where p(t) = —dP{t)/dt. Before applying (9), we must also ensure that 
a and p are small so that each visit to the accelerator mode is uncorrelated with the previous 
one. (If (T or p is not small, the particle may be immediately scattered back into the trapped 
region after being scattered out of it, instead of being swept far away from the island. This 
effectively increases the area of the trapped region. Such considerations were not necessary 
for the large-angle scattering model.) Taking p = 0, we have p{t) ^ q(at/A'^) where q is 
some function which, for a given K, applies for all first-order accelerator modes and which 
is independent of a. A is the scale length of the island which from (12) is proportional to 
1/n. Here again, a must be small to be able to write p{t) in this way because we need to 
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bc able to separate the slow time scale of the motion due to the noise in (2) from the fast 
time scale due to the standard mapping itself. This form for p(t) applies for any diffusive 
noise model if we regard u as a measure of the extent of the Green's function response after 
a unit of time under the action of the noise alone. Using (9), we find that, when e satisfies 
(13), Dig ~ l/((7e^) and Dig/Dqi ~ l/iae^). This is confirmed by fig. 2, where we see that 
D/Dqi decreases by roughly a factor of 16 when e is doubled. 

5. Conclusions 

When noise is added to the standard mapping the diffusion coefficient consists of two 
parts. One part is primarily due to short-term correlations in the stochastic region of phase 
space and this part exhibits nearly sinusoidal oscillations [2] . The other part is proportional 
to the inverse of the noise parameter a and exists only when accelerator modes are prescmt. 
When the noise is weak, this latter part dominates. The role of the noise is to ensure that 
every trajectory eventually visits the accelerator modes. This effect has been confirmed 
numerically for the first-order accelerator modes. These modes exist in windows given by 
(13) and their relative importance decays as for diffusive noise (2) and as for large- 
angle scattering (3). Although many higher-order accelerator modes exist (see table I), they 
are so small that they do not contribute significantly to D for the values of a that can be 
dealt with computationally. In the absence of noise, the diffusion coefficient is not fully 
defined until the ensemble in (4) is specified. However, anomalies in the diffusion coefficient 
are seen even when the ensemble in (4) is taken to be trajectories in the stochastic region. 
This is apparently because trajectories can be greatly accelerated while in the vicinity of 
the accelerator modes. 

We showed in section 4 that the accelerator modes normally exhibit a universal behav- 
ior. Depending on the noise model, the diffusion coefficient therefore exhibits some universal 
features as the parameter e is varied. 

An interesting accelerator mode which is not accurately described by (11) is found 
near e = 3. At this value of the parameter, Chirikov [1] noticed that the distribution of 
particles after 100 iterations of the standard mapping was characterized by two diffusion 
rates. (See fig. 5.8 of ref. 1.) This is due to the existence of a 4th order accelerator mode 
with acceleration \a\ /{2tt) = j. The history of this mode is somewhat unusual. The mode 
is born as a tangent bifurcation at e = 2.9453. The period-4 orbit becomes unstable at 
e = 2.9775, giving rise to a pair of period-8 orbits. However, these period-8 orbits do 
not become unstable as e is increased. Rather they are reabsorbed into the period-4 orbit 
making it stable again (this happens at e = 2.9950). At e = 3.1068, it becomes unstable as 
an ordinary hyperbolic orbit. Simultaneously, two stable period-4 orbits are created. They 
become unstable at e = 3.1449 giving rise to twice as many period-8 orbits. The standard 
period-doubling process now takes place as e is increased. (So actually the accelerator mode 
exists for e somewhat larger than 3.1449.) This degenerate behavior is probably due to the 
fact that this accelerator mode lies on one of the symmetry lines of the standard mapping, 
V = nn + iesin(a;) with n being an integer. 

Similar effects are expected in other mappings which allow accelerator modes. Such 
mappings are ones which are periodic in the velocity direction. For example, an anomalously 
large diffusion coefficient was found at a certain value of the parameter for the mapping stud- 
ied in ref. 7. (The divergent behavior of the difl[usion coefficient for this case has been noted 
by Antonsen and Ott [9].) However, such mappings are often derived as approximations to 
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mappiiigs which arc not pc;riodic in the velocity direction. For instance (1) may be derived 
from similar mappings in which e is effectively a function of v, e{v) [14]. These more gen- 
eral mappings do not have accelerator modes. Nevertheless, if the dependence of e on t; is 
weak, we still expect there to be some effect due to the presence of the "latent" accelerator 
modes. These are channels in phase space which allow particles to be uniformly accelerated 
from a velocity vq such that e{vo) w eo to a velocity V2 such that €{v2) ~ €2 (eo and €2 
being the parameter values for the birth and death of an accelerator mode in the standard 
mapping). No noise is required to bring particles into these channels. This short cut will 
result in an enhancement of the diffusion coefficient between vo and V2- Diffusion on longer 
velocity-space scales will still be limited by the slower diffusion elsewhere. 

Finally, we note that the 1/cr dependence of D may arise in the study of particle 
confinement in fusion devices, whenever a small group of particles is unconfined in the 
coUisionless limit. (For problems involving diffusion in real space, the unconfined particles, 
for which the position coordinate is monotonically increasing, correspond to accelerator 
modes.) Collisions, which can scatter these particles into confined orbits, serve to slow the 
loss of these particles leading to a diffusion process which is inversely proportional to the 
collision frequency [15]. 
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Tables 

Table I. First- and second-order stable accelerating fixed points for the standard mapping 
(1). Here, N is the order of the fixed point and a is its mean acceleration. The subscripts 
and 1 denote the values of e, a;, and v when the fixed point first appears and when it becomes 
unstable. These fixed points together with those generated from these by the symmetries 
of (1) give all the first- and second-order accelerating fixed points for |e| < 20. The entries 
labelled with asterisks result when the preceding first-order point becomes unstable. 
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Fig. 1. The dependence of D/Dgi on a where Dqi = ■ The bars indicate the error (± 
one standard deviation). The lower two curves are for the diffusive noise model (2) with 
p = and (a) e = 6.6 and (b) e = 12.8. The upper curve (c) is for the large-angle scattering 
model (3) and e = 6.6. The dashed lines show the l/cr dependence. The point with the 
lowest value of a on each curve was calculated with fcT = 2 x lO'^ and kL = 10^. Lower 
values of these parameters were used for the larger values of a. The point on the right (d) 
gives D/Dqi for cr = and initial conditions in the stochastic region. Here kT = 10® and 
fci = 5 X 10^ 
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Fig. 2. D/Dqi as a function of e for the standard mapping with diffusive noise (2) with 
p = and (T = 3 X 10~^. The bars show the error as in fig. 1. Here kT = 6 x 10^ and 
fci = 3 X 10'*. The dips at e « 6.65 and 12.75 are due to the fourth-order resonance while 
those at e « 6.95 and 12.9 are due to the third-order resonance. 
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Fig. 3. Comparison of the quadratic mapping (11) with the standard mapping (1) in the 
regions of accelerator modes. The first figure (a) shows trajectories for the mapping (11) 
with K = 0.265 with 9 different initial conditions. There is a first-order island at the origin 
whose center is /i. The other two largo islands are part of a chain of 4 fourth-order islands. 
The other figures (b and c) show trajectories for the mapping (1) with e = 18.9620 and 
15.0425 where first- and second-order accelerator modes exist. The coordinates {x, v) in 
(1) have boon transformed to {X^Y) according to the proscription given in Section 4 with 
K = 0.265. I.e., for (b) the transformation is (12) with n = 3; and for (c) it is x/{2tt) = 
0.2254 - 10-3(3.4298X -|- 0.0019y -|- O.OOOIA:), v/{2tt) = 0.1826 - - 10-3(3.9784X - 
0.272iy - O.Oligii'), and e = 15.0416 -|- 3.6453 x lO'^K. 



-18- 




,L. I 1- - 1 i 1 

'O 0.25 0.50 0.75 1.00 1.25 




Fig. 4. The trapped region of phase space for the quadratic mapping (11). Plots (a) and 

(b) shows the intersection of the trapped region with the two symmetry hnes V = K — X"^ 
and V = Q. The trapped region is represented in black. These plots are calculated with a 
resolution of 0.002 in X and K. Trajectories of length 10 000 were used in the test for being 
trapped. The curve (c) shows the total area A of the trapped region. This was calculated 
with a resolution of 0.005 in X and Y and at intervals of 0.01 in K. Here the trajectory 
length was 1000. 



